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ABSTRACT 



We use the high magnification event seen in the 1999 OGLE campaign light curve of image C of the quadruply imaged gravitational 
lens Q2237+0305 to study the structure of the quasar engine. We have obtained g'- and r'-band photometry at the Apache Point 
Observatory 3.5m telescope where we find that the event has a smaller amplitude in the r'-band than in the g'- and OGLE V-bands. 
By comparing the light curves with microlensing simulations we obtain constraints on the sizes of the quasar regions contributing to 
the g'- and r'-band flux. Assuming that most of the surface mass density in the central kiloparsec of the lensing galaxy is due to stars 
and by modeling the source with a Gaussian profile, we obtain for the Gaussian width 1.20 x 10" s/M/OAMq cm < CTgi < 7.96 x 
10'^ VM/O.lMo cm, where M is the mean microlensing mass, and a ratio cr,-i jcTgi = 1.25^q jj. With the limits on the velocity of the 
lensing galaxy from Gil-Merino et al. (2005) as our only prior, we obtain 0.60x 10'' ^IMIQAMq cm < cTj' < 1 .57x 10'^ ^IMIQAMq cm 
and a ratio a^' IcTgi = 1.45^j5 '5 (all values at 68 percent confidence). Additionally, from our microlensing simulations we find that, 
during the chromatic microlensing event observed, the continuum emitting region of the quasar crossed a caustic at > 72 percent 
confidence. 
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1. Introduction 

■ One of the best studied quasar lensing systems yet is 
[ Q2237 -1-0305. This quasar was disc overed by the Cent er for 
. Astrophysics (CfA) redshift survey dHuchra et alj Il985l) . The 
system is made up of a barred spiral galaxy (z=0.0394) and four 
images of a quasar at a redshift of z- 1 .695 that appear in a nearly 
symmetric configuration around the core of the spiral galaxy at 
a distance of approximately 1" from the center. 

Q2237 -1-0305 is an ideal system to study quasar microlens- 
' ing. At the position of the quasar image s the optical depth 
' due to micro lensing by the stars is high ( Kavser et a l.l [19861 
iKavser & Ref sdal 1989i: IWambsganss et al.lll990l) . Furthermore 
the expected time-delay betw een the images of this quasar is 
of the order of a day or less ( Schneider et al.l Il988t iRix et al] 
119921: IWambsgans s & Paczvnskilll99i r As microlensing events 
for this system are of the order of months or years, it is easy to 
distinguish intrinsic variability from microlensing. Another big 
advantage is the fact that the quasar is ~10 times farther from 
us than the lensing galaxy. This leads to a large projected veloc- 
ity of the stars in the source plane and thus a short t imescale of 
caustic events. Not long after its discovery in 1985. llrwin et al.l 
(119891) observed the first microlensing signature in this quasar. 

The quasar has been monitored for alm ost two decades 
by different surveys (e.g., Corriganet al. 1991: |Pe"n et al.|[T993t 



lOstensen etai]ll996l: I Alcalde et all 120021: LSchmidt et al.ll2002h . 
The one that has the longest and best sampled data is that 
of the Optical Gravita t ional Lensing Exper iment (OGLE) team 
dWozniak et all l2000t lUdalski et all |2006). They have moni- 
tored Q2237H-0305 since 1997 delivering the most complete 
light curves for this system in which many microlensing events 
can be seen. These data ha ve been used f or various studies 
on the system: I'Wvithe et al.1 Jlo OOa b.d e. j): lYoneharal ( 1200 lb : 
Shalvapi n et alj ( |2002 |) and 'Koc hanekl ([2004) used this data to- 
gether with microlensing simulations not only to obtain limits on 
the properties of the quasar such as size and transversal velocity, 
but also on the mass of the microlensing objects. 

The OGLE team has monitored this quasar with a single V 
band filter. However, IWambsganss & Paczvnskil (1199 lb showed 
that color variations should be seen in quasar microlensing 
events and how they would provide additional information on 
the background source. We know from thermal accretion disk 
models that emission in the central r egions of a quasar should b e 
bluer than in the outer regions (e.g.. lShakura & Sunvaevl[T973h . 
As an approximation we can assume the central engine of the 
quasar to be of the order of 1000 Schwarzschild radii: 
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* Based on observations obtained with the Apache Point Observatory 
3.5-meter telescope, which is owned and operated by the Astrophysical 
Research Consortium. 
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where r, is the Schwarzschild radius, R.^^ is the radius of the 
accretion disk and Mbh is the mass of the black hole. We can 
compare this value with the Einstein Radius (Re) - the length 
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Fig. 1. APO g' image of the quasar system Q2237+0305. 

scale in q uasar microlensing - w hich can be given in the source 
plane by jSchneider et alJll992h : 

where we have replaced the values of the distances D/, Dg and 
Dis to the ones corresponding to Q2237+0305. The accretion 
disk scales are thus comparable to those of the Einstein rings of 
stars in the lensing galaxy. This shows that differential magni- 
fication of the different emission regions of the quasar can be 
observed. 

Some mu lti-band observation s of this system hav e been 
carried out (Corrigan etaD Il99ll: IVakuhk et all 1 19971 l2004t 
[Koptelo va et al. 2005), however most of them have either large 
gaps between the observations, or no obvious microlensing 
event. 

A very int eresting event, that has been o ne of the bases 
for prediction s (Wvithe et al. 2000b) and study (Yonehara 2001; 
IShalvapin eTa l. 2002; Vakuhk et al. 2004), is the high magnifi- 
cation event of image C in the year 1999 as observed by OGLE. 

In this paper we present data of this particular event obtained 
in two filters at Apache Point Observatory (APO). We analyze 
the event using both the OGLE light curve and the APO data 
and study the implications for the size of the quasar emission 
region using microlensing simulations. We use a flat cosmology 
with Q„, = 0.3 and Hq = 70 km s'^Mpc^K 

2. Apache Point observatory Data, Observations & 
Data reduction 

For this study we use two different datasets. The first one is that 
of th e OGLlJ] team light curves for this quasar dWozniak et alj 
1200 0; Udalski et al. 2006) for which we select the data between 
Julian days 2451290 (April 21, 1999) and 2451539 (December 
26, 1999), which comprises 83 data points. The result is a dense 
light curve for the event we analyze. The second data set is 
APO two-band monitoring data for Q2237 -1-0305. We selected 
the nights between May 17 1999 (JuUan day: 24511316) and 
January 8 2000 (Julian day: 24511551) which also track the 
event, comprising 8 nights with data points. 

' http://bulge.princeton.edu/~ogle/ogle2/huchra.html 



The APO data were taken with the 3.5 m telescope at Apache 
Point Observatory using the Seaver Prototype Imaging camera 
(SPIcam) in both the SDSS g' and r' filters. The SPIcam has 
a 2048x2048 pixels CCD and a minimum pixel scale of 0.141 
"/pixel, however, for the observations we use, the pixels were 
binned to 0.282 "/pixel. 

2.1. Standard CCD Reduction 

For the standard CCD reduction we use a combination of the 
astronomical image reduction software package iraf and idl data 
language. We apply bias correction, create one master flat-field 
image for each night and each filter and calibrate the science 
frames with these. Cosmic rays are eliminated from the science 
frames using median filtering outside the sources. 

In order to remove the remaining cosmic rays and imperfec- 
tions of the CCD as well as to improve the signal to noise ratio 
of the images, we combine the frames corresponding to the same 
night, excluding those with asymmetric point spread functions 
(PSFs). In the end we obtain a clean high quality image for each 
night in g' and r' filters. 

2.2. Photometry 

The photometry for this system is made difficult because 
the quasar images are seen through the core of the lens- 
ing galaxy. Our first attempt, following the procedure as in 
Wozniak et al. (2000), was to use an image subtraction method 
( Alard & Luptoni .1998.) . However our field of view lacks the 
high number of bright stars required in order to run a suc- 
cessful correlation and PSF character ization for the procedure. 
This is why we use galfit (Peng et al. 2002) instead, galfit is a 
galaxy/point source fitting algorithm that fits 2-D parameterized 
image components directly to the images. 

2.2.1 . GALFIT model 

For the galaxy we use the m odel shown bv lSchmi"d3 ( Il996h (see 
also iTrott & Websted |2002|) from HST F7815LP band data, in 
which the lensing galaxy is parameterized with a de Vaucouleurs 
profile for the bulge and an exponential profile for the disk. The 
bulge's de Vaucouleurs profile is constrained to a 4.1" (3.1 kpc) 
scale length, a 0.31 ellipticity and a 77° position angle. For the 
exponential profile we use a 1 1.3" (8.6 kpc) scale length, a 0.5 el- 
lipticity and a 77° position angle. These profiles are constrained 
to have the exact same central position within the image and a 
fixed magnitude difference. 

The quasar images are parameterized as point sources. The 
relative separation bet ween them is fixed in the galfit code to 
the values obtained by Blanton et alj ( Il998l) from UV data and 
the separation between the group of images and the center of 
the bulge is fixed to the values obtained with HST F7815LP fil- 
ter data (HST proposal ID 3799): Ara = -0.075" and Adec = 
0.937". 

The PSF we use to convolve the different parameterizations 
is taken from the bright star (STD) shown in Fig. [T] After fit- 
ting the different profiles, the brightness of each component is 
measured. We run the galfit routine over each night in each one 
of the two bands available. The residuals from this fitting tech- 
nique are low, usually of the order of 5 percent. In order to ob- 
tain error bars on the magnitude measurements, we create 500 
Monte-Carlo realizations of the observed images and determine 
the measurement uncertainty from the scatter of the galfit re- 
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Fig. 2. Top panel: Black squares show the APO g'-band light 
curve (relative to the faintest data point) of the high magnifica- 
tion event seen in image C of Q2237+0305 in 1999. For compar- 
ison, in light gray, the OGLE V-band light curve sampling this 
event is shown. Middle panel: APO r'-band light curve for the 
event (black squares) and the OGLE V-band light curve plotted 
in light gray. Bottom panel: APO g'-r' color curve. The dashed 
line is shown to guide the eye. 



suits. The g'- and r'-band light curves for image C are shown 
in Fig. [21 The measurements for images A and C are given in 
Table [T] No reliable photometry is obtained for images B and D 
because they were too faint. 

As shown in Fig. |2l the APO g'-band light curve shows a 
very similar behaviour to the OGLE V-band light curve, how- 
ever, the APO r'-band light curve shows a lower amplitude at 
the brightness peak. Thus, the resultant g'-r' color curve, does 
not appear flat but shows a chromatic variation during the mi- 
crolensing event. 



3. Microlensing Simulations 

In order to compare the observed data with simulations we gen- 
erate source-plane magnification patterns fo r quasar images C 
and B using the inverse ray shooting method (Wambsgan ss et all 
[1990; Wambsganss 1999). We assume identical masses for 
all the microlenses distributed in the lens plane (the results 
do not depend on the mass distribution of the microlenses 



Table 2. Local macro-lensing p arameters for quasar images B 
and C from lSchmidt et aD(ll998l) . k is the scaled surface density, 
y is the shear, fitot is the total magnification. 
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0.36 


0.42 


4.29 
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0.69 


0.71 


2.45 



(iLewis & Irwin|[T996h ). Approximately 10^' rays are shot, de- 
flected in the lens plane and collected in a 10,000 by 10,000 
pixels (equivalent to 50 Einstein radii by 50 Einstein radii) array 
in the source plan e. The values of sur face density k and shear 
y are taken from Schmidt et all ([1998 ) (Table [2]). The conver- 
gence K is separated into a compact matter distribution and a 
smooth matter distribution Kc- Even though we explore different 
stellar fractions, as the quasar images are located within the cen- 
tral kiloparsec of the lensing galaxy, we focus our work mainly 
on the — = 1 stellar fraction. 

As shown in Sect. 1, in quasar microlensing the source size 
is an important parameter b ecause it i s comparable with typical 
caustic scales (e.g. Kavser et al.lll986l) . In order to study a sam- 
ple of different source sizes, th e raw magnification p atterns are 
convolved with a set of profiles- TMortonson et al.l (12005) showed 
that microlensing fluctuations are relatively insensitive to the 
source shape, so it is of little conseq uence whether we ch oose 
a standard model accretion disk ( Shakura & Sunvaevll 19731) or a 
Gaussian profile. For simplicity we choose a Gaussian profile for 
the surface brightness where the extent of the source is described 
by the variance cr. The full width half maximum (FWHM) is de- 
fined as 2.35cr. 

We vary the FWHM of a Gaussian profile from 2 to 120 pix- 
els (which corresponds to 0.01 Er and 0.60 Er, respectively) for 
quasar image B and from 2 to 216 pixels (which corresponds to 
0.01 Er and 1.08 respectively) for quasar image C (the addi- 
tional magnification patterns for image C are used for the color 
curve fitting explained in Section 5). For patterns below 12 pix- 
els the step size is 2 pixels, and for patterns above 12 pixels the 
step size is 4 pixels. We find that this linear sampling gives more 
accurate results when interpolating compared to a logarithmic 
sampling. 

For a specific convolved pattern we can now extract light 
curves. By defining a track (path of the source in the source plane 
magnification pattern) with a starting point, direction and veloc- 
ity, we extract the pixel counts in the positions of the pattern 
defined by this track using bi-linear interpolation. The values 
are then scaled with the magnification values of Table[2lfor each 
image and converted into magnitudes. 

4. Light curve fitting 

4.1. OGLE V light curve 

Since quasars vary intrinsically, we need to separate this vari- 
ation from microlensing. We therefore study the difference be- 
tween two light curves extracted from microlensing patterns for 
different images to eliminate this effect and allow for an ad- 
ditional magnitude difference (the time delay is negligible, see 
Sect. 1). 

Using the V band OGLE light curves sampling this event 
in both images C and B of the quasar ( Wozniak et al. 200(| 
Udalski et al. 2006) we construct the difference light curve C-B. 
Simulated light curves for each image are extracted using tracks 
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Table 1. APO light curves for quasar images A and C. Magnitudes are shown relative to the value of the last data point of image C 
in the g' filter. 



Julian Day 
+2451000 


A [mag] 


g' 

err [mag 


A 

A [mag] 


r' 

err [mag] 


A [mag] 


g' 

eiT [mag] 


C 

A [mag] 


r' 

err [mag] 


FWHM 

[arcsec] 


1315 


-0.822 


0.016 


-0.675 


0.012 


-0.533 


0.016 


-0.381 


0.013 


1.6 


1334 


-0.826 


0.056 


-0.680 


0.021 


-0.554 


0.056 


-0.393 


0.021 


1.4 


1373 


-0.925 


0.017 


-0.729 


0.013 


-0.672 


0.017 


-0.445 


0.013 


1.2 


1385 


-0.971 


0.017 


-0.782 


0.015 


-0.635 


0.017 


-0.411 


0.015 


1.0 


1400 


-0.912 


0.019 


-0.757 


0.014 


-0.487 


0.019 


-0.298 


0.014 


1.1 


1531 


-1.148 


0.024 


-0.887 


0.015 


-0.013 


0.026 


0.067 


0.015 


1.1 


1541 


-1.158 


0.015 


-0.909 


0.015 


-0.037 


0.016 


0.016 


0.018 


1.1 


1551 


-1.137 


0.018 


-0.916 


0.014 


0.000 


0.024 


0.094 


0.015 


1.2 



over microlensing patterns created independently for each im- 
age and convolved with a Gaussian profile with the same size. 
By repeated light curve extraction with variation of the parame- 
ters, the best fit to these are obtained from the comparison with 
the OGLE C-B difference light curve. What is actually fitted are 
the parameters that define the track from which both light curves 
are obtained. This track is constrained to have identical direction 
and velocity in both the patterns C and B, but the starting point 
can be different among the different patterns. It is important to 
remark that the direction is set to be the same in both patterns, 
taking into consideration that the shear direction between images 
C and D is approximately perpendicular ( Witt & Mao 1994). 

For the minimization method we use a Levenberg-Marquardt 
least squares routine in idi0. Th is routine is a ;t^^-based minimiza- 
tion based upon MINPACK-1 (iMore et al.1 1980). As we are us- 
ing 68 percent values for our errors, the is simply calculated 
as: 




where DOF is the number of degrees of freedom, y is the mea- 
sured value, X is the independent variable model and (7^" is the 
one sigma error of the particular measurement. In the case in 
which we fit the OGLE data for this event, the number of degrees 
of freedom is 75. y-, are the differences between the OGLE light 
curves for images C and B, cr"'' are the uncertainties measured 
by the OGLE team, .jc, are Julian days and /(x,) are the differ- 
ence between of the two light curves extracted from the source 
plane magnification patterns for images C and B convolved with 
an identical source profile. We also include a parameter which 
we call mo that accounts for the magnitude offset between the 
images. 

As the size of each one of the convolved patterns is very 
large compared to the time scale of the chosen event (of the or- 
der of several hundred times), and considering the fact that it is 
relatively easy to find good fits because of the huge parameter 
space that such large patterns give, we need to determine a high 
number of tracks for each pattern. To obtain a statistical sam- 
ple, we search for 10,000 tracks for each considered source size. 
The first guesses for each one of the tracks are distributed uni- 
formly. However, during the fitting process the starting points of 
the tracks are constrained to stay within a square of 200 pixels 
sidelength around the random initial value in order to force a 
sampling of the whole pattern. After this minimization is done 
we have a track library containing 3x10^ fitted tracks (10,000 

^ http://cow.physics.wisc.edu/~craigm/idl/ 



tracks for 30 different source sizes with best-fitting velocity and 
magnitude offset). 

4.2. APO Color Curve 

Fig. |2] shows that the OGLE V-band and the APO g'-band light 
curves are very similar. Since the microlensing effect depends 
strongly on the source size, we assume, for simplicity, that the 
OGLE V-band continuum region is of similar size as that of the 
APO's g'-band. Here we determine the ratio of the APO r'-band 
region and the g'-/OGLE V-band regions that follows from the 
color curve in Fig.|2l 

While the parameters that define the tracks are fixed by the 
procedure shown in the previous section, the two parameters that 
influence the shape of the color curve in this step are the ratio 
between the sizes of the two regions and the intrinsic color g'-r' 
of the quasar (cq). These are the only parameters allowed to vary 
during this step. As before, we use the patterns convolved with 
Gaussian profiles up to 120 pixels FWHM for image C to obtain 
the g'-band light curve. To obtain the r'-band light curve we use 
the full set of patterns for image C. During the fitting process we 
interpolate the magnitude values of the extracted light curves in 
order to obtain continuous values for the source size ratio. After 
this second fitting procedure, the values of the track library 
are updated: ^2 =^2^^^ + 

5. Considerations 

5.1. Degeneracies 

As source size and microlens mass are degenerate, we express 
the size values in Einstein radii or scaled by a stellar mass as 
shown in equation (|2]i. Another important degeneracy is the one 
between the size of the source and the transversal velocity we 
fit. If we fix the path velocity of a source and increase the size of 
it, its light curve will get broader Conversely, if we fix a source 
size and decrease the velocity, the light curve will also become 
broader because it takes longer time for the source to cross the 
magnification region. 

5.2. Velocity considerations 

The velocity one measures for a microlensing event is composed 
of velocities of the observer, the lenses and the source. The 
Earth's motion relative to the microwave background is likely 
l :o be almost p arallel to the direction of the quasar Q2237 -1-0305 
dWitt & Mad 1 19941) so the velocity of the observer is negligi- 
ble compared to the other velocity values. We assume the phys- 
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Fig. 3. Example track for a 0.05 Er (10 pixel) FWHM g'-band 
source size and a source size ratio (Tr'/o-g' = 1.48. In the top 
panels the fitted track is shown in a 1 .QEr x 1 .QEr section of 
the microlensing patterns. In the middle panel we plot the fitted 
difference light curve (solid line) and the observed OGLE differ- 
ence light curve. In the lower panel the best-fitting color curve 
(solid line) and the APO g'-r' color curve are shown. 

ical velocities of the lens and the source of the same order of 
magnitude, thus, the main contribution to the total velocity is 
the velocity of the lens due to the fact that the source veloc- 
ity is we ighted by ~3.5 percent compared to the lens velocity 
jKavser e t al. 1986). 

The velocity of the lens has two components: one of 
the bulk velocity of the galaxy and the velocity disper- 
sion of the stars in the galaxy (microlenses) v^. Since the 
time span of the event that we fit is ~200 days, the sig- 
nificance of individual st ellar motions in the fi tt ing sh ould 
be fairly low. Moreov er, Kundic & WambsganssI (1 1993b and 
IWambsganss & Kundid d 1995 1) show that, in a statistical sense, 
the overall effect of the individual stellar motions can be approx- 
imated by an artificially increased velocity (a factor » 1 .3) of the 
source across the pattern. 

The scale we use for fitting the velocity is pixels per Julian 
day. Using equation (|2]i we can see that this turns into: 



pix 
Id 



0.005 



jd 



32870 - 



I M km 
(O.lMoT 



(4) 



in the source plane. Using the appropriate velocity scaling 
dKayseret aL,1986) . we find: 1 ^ ~ 3000 ^ in the lens plane. 



5.3. Size Considerations 

In a standard thin accretion disk model of a quasar the black 
hole is surrounded by a thermally radiating accretion disk 



(IShakura & SunvaevI [1971 . iKochanek et"!!] (l2006h show that 
the radii where radiation with wavelength A is emitted scale 
as /l3 . Thus: 



(5) 



For the observed data from APO, the transmission curve for the 
SPIcam g' and r' filters have central wavelengths of ~4700A 
and ~6400A, respectively. Using these values and equation (|5]) 
we obtain a predicted thin-disk ratio between the r' - and g'-band 
source sizes of ~ 1.5. 



6. Statistical Analysis 

There are different ways to approach quasar microlensing 
studies. It ca n be by analytical model fitting of the light 
curves (e.g., Yonehara' '2001), studying the structure func- 
tion (e.g., Lewis & Irwin 1996), doi ng statistics of parame- 
ter v ariations over time intervals (e.g., Schmidt & Wambsganss] 
11998: Wambs ganss et al. 2000: Gil-Merino et al., 2005) . obtain- 
ing probability dist ributions with hght curve derivatives (e.g 
IWvithe et al.'2000c') orBayesian analysis (e. g..l Kochanekl2004 
Similar to Kochanek (2004), we construct a probability distribu- 
tion for the quantities of interest (V/g'-band source size, source 
size ratio and transversal velocity) from our track library. 

Each one of our tracks is a fit to the light and color curves 
and therefore has a value attached to it. We therefore infer 
information from the whole ensemble of models: the track li- 
brary. Using the standard approach for ensemble analysis (e.g.. 



iing 

ISambridgdl 19991) . we assign a likelihood estimator to each track 

f,-: p(ti) oc exp(^Y^). Using this likelihood estimator the statisti- 
cal weight for each track can be written as: 

P(ti) 



w(ti) = 



n{ti) 



(6) 



where w(f,), p(tj) and n(f,) are, respectively, the weight, the like- 
lihood obtained from the and the density parameter for each 
track fj. The density parameter describes the local track density 
in a multidimensional grid in which each dimension corresponds 
to a parameter of interest. The density is normalized so that the 
sum of all the density parameters in the grid equals one. By sum- 
ming over the statistical weights, we can calculate probability 
distribution histograms for all the parameters of interest. 

Gil-Merino et aL (.2005) (from now on G-M05) found an 
upper Umit of 625 ^ (90 percent confidence) on the effective 
transverse velocity of the lensing galaxy in Q2237-I-0305. They 
performed microlensing simulations assuming microlenses with 
M = 0.1 for three different source sizes yielding similar re- 
sults. We use their probability distribution for the velocity ob- 
tained with the largest source size value (Fig. 5 in G-M05) in a 
further analysis step, where we factor it into our own probability 
distributions by importance sampling. In other words, we scale 
the probability of finding a track based on the G-M05 prior, and 
then reobtain the best-fitting value for each parameter and confi- 
dence regions. 

7. Results and Discussion 

As described in the previous section, we obtain limits on the 
source size and transverse velocity from two sets of probability 
distributions. The first corresponds to the parameters obtained 
without using any prior, and a second one in which we apply 
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Table 3. Results and confidence limits for g'-band source size, ratio between r'- and g'-band source size and transverse velocity 
(projected to the lens plane). 

No Prior With G-M05 Prior 

Source Size o-g> Ratio Velocity (*) Source Size (Tgi Ratio Velocity 



M 
O.lMo 



xlO'^[cm] [km/s] xlO'-^[cm] [Icm/s] 

SAO^lfg 4.60^_ll^ l-25!!Jt^ 2930 2.34+»« l-33!»-^f| 1.45:«* 682+^27 



(*) The velocity obtained without the G-M05 prior is shown without error bars as we do not obtain limits for it. 



the G-M05 prior on the velocity of the lens galaxy. For each 
probability distribution of a particular parameter we select the 
68 percent confidence levels. This is done by making horizon- 
tal cuts starting from the highest probability and screening down 
until 68 percent of the cumulative probability is reached. The 
distribution, the best-fitting values (or the average of the 68 per- 
cent confidence region where no single best-fitting solution is 
found) and the 68 percent confidence limits for the relevant pa- 
rameters are shown in Fig. |4] and Table [3] After applying the 
G-M05 prior only ~1000 tracks carry 99 percent of the statisti- 
cal weight. This makes the probability histogram for the source 
ratio sparsely populated for cr,.>/crg' > 2. Therefore, we require 
at least 70 of these 1000 tracks to be part of each bin in the his- 
togram. 

Without any velocity prior (see Fig. HJi, we have more fast 
tracks than slow tracks. No limit on the transversal velocity 
can be determined. This mostly is due to the degeneracy be- 
tween source size and transverse velocity (see Section 5.1). 
Using the G-M05 prior, we obtain a best-fitting velocity of 

682^3^g -^J o fif^ k-tn/ s (projected into the lens plane). Note, how- 
ever, that the distribution shows a longer tail towards the higher 
velocities when compared to the G-M05 probability distribution 
because the steepness of the event fitted prefers fast tracks or 
small source sizes. 

The average Gaussian width we obtain on the OGLE V/APO 
g' source size without setting the velocity prior is ctj/ - 

S.lO+lfg X 10-^ Re or 4mtll^yA0\f^cm. These limits 

agree with those obtained by Yoneharal (|2001|) and lKochanel3 
(I2OO4I) (who used a logarithmic prior on the transversal ve- 
locity). When we impose the velocity constraints of G-M05 
we obtain an OGLE V/APO g' source size with a Gaussian 
width of o-g, - 2. 34+ J 28 ^ ^O'^Re which is equivalent to 

1 .33;^g73 X 10'^ sj o.ute '^^- '^^^^ makes the upper limit tighter by 
a factor of five compared to the previous result, pla cing it just 
below the 68 percent confidence limit obtained by iKochaneM 
(Hool, very close to their resolution limit. For niQ we obtain 
-1.75 < mo < 0.25 without and -0.75 < mo < 1.25 with the 
G-M05 prior, respectively (see Section 4.1). 

The source size ratio between the r'- and g'-band emit- 
ting regions of the quasar is cTr'/crg' = 1 .25'^^^''\^^ without, and 
CTr'/cTg' = 1.45;^Q 25 ^ith the G-M05 prior, respectively, coupled 
with an intrinsic color parameter: cq = -0.18 ± 0.05 in both 
cases (see Section 4.2). As shown in Fig. HI both source size 
ratio distributions are similar. The measurements are close to 
the theoretical value (f~1.5) estimated in Section 5.3 from the 
[Shakura & Sunyaev ( 1973) thermal profile. 

For a long time there has been disc ussion on the nature 
of the image C event in year 1999 (e.g.. IWyithe et allllOOObl 



Table 4. Probabilites of the image C high magnification event 
being produced by the source crossing a caustic and of the g' 
and r'-band region touching a caustic. The table shows the ob- 
tained values for the pr obability both using and not using the 
iGil-Merino et alj (l2005b velocity prior. 





No Prior 


With G-M05 Prior 


cross 




cross g' r' 


0.75 


0.92 0.97 


0.72 0.89 0.94 



lYoneharallMTt IShalvapin et aLll200l iKochanekl 12001 . Using 
different methods it was found that this event is likely to have 
been produced by the source passing near a cusp. We investi- 
gate this for the tracks in our library by computing the location 
of the caustics for each of the magnification patterns using the 
analytical method by Witt (1990, 1991) (combination of both 
ray shooting simulation s and analytical caustics are shown in 
IWambsganss et al.lll992l) . 

For each track we determine whether either the center, the 
centericTg' or the center±cr,./ of the source touches a caustic at 
any time. Our analysis agrees with the previous results when we 
use the OGLE light curve to fit for this event. However, to re- 
produce our APO color dataset, tracks that cross a caustic are 
favored (see Table |4]l. We find that a source with a radius as 
big as cTr' touched a caustic with 97 percent and 94 percent 
of confidence without and with the G-M05 prior, respectively. 
Furthermore, the source, regardless of dimensions, crossed a 
caustic with a 75 percent and 72 percent of confidence with- 
out and with the G-M05 prior, respectively (i.e. the center of the 
source touched a caustic). 

The probabilities of the source touching a caustic at different 
radii are bigger in the case in which we do not use the G-M05 
prior as seen in detail in Table |4] This is due to the fact that a 
higher velocity makes a track cover more distance through the 
pattern and thus is more Ukely to encounter a caustic. 

8. Conclusions 

We present two band (g' and r') APO data covering the high 
magnification event seen in image C of the quadruple quasar 
Q2237-I-0305 in the year 1999. We find that the amphtude 
of the brightness peak is more pronounced in the g'-band 
than in t he r'-band. This is also consistent with the observa- 
tions by IVakulik et al.l (l2004t) (see also the recent analysis by 
lKoptelovaetal.ll2007t) . 

By using this data together with the well known OGLE data 
( Wozniak et al. 2000; Udalski et al. 2006) and combining it with 
microlensing simulations, we have been able to obtain limits on 
the size of regions of the quasar's central engine emitting in these 
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1.45+°25. Both values for the ratio be- 



bands: Gaussian width o-g- = 4.60+^-;^J;x 10'^ yfMjTlM^ cm and 

Because of the degeneracy between source size and trans- 
verse velocit y we use a prior on th e velocity obtained from 
the work of ! Gil-Merino et alj ( l2005l) in order to improve the 
results: Gaussian width cr, - ^ ^^+024 
cm and cr^'/o-g/ 

tween the source sizes are c l ose to the ratio obtained for 
a face-on IShakura & SunvaevI (1973h accretion disk (f~1.5) . 
Recent s tudies (e.g. JPoindexter et alj|2007 ': Morgan et al. 2007^, 
iPoolev et al. 2007) suggest microlensing yields a slightly bigger 
value for the ratio than that obtained analytically with the thin 
disk model (see Section 5.2), also in agreement with our results. 

We also show that this event was probably produced by the 
source directly interacting with a caustic as we obtain probabili- 
ties of 97 percent and 94 percent that the r'-band emitting region 
touches a caustic without and with the G-M05 prior, respectively, 
and of 75 percent and 72 percent that the source center (regard- 
less of size) crossed a caustic without and with the G-M05 prior, 
respectively. 
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Fig. 4. Probability histograms for g' or V source size, r'/g' source size ratio and transverse velocity (projected to the lens plane), 
respectively. The dotted Une shows the best-fitting values (or average of the 68 percent confidence region if no single best-fitting 
value is found). The error bars describe the 68 percent confidence limits. 



